double erf(double x);